Wobbling excitations and tilted rotation in Lu 
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Using a microscopic self-consistent model, we analyse wobbling excitations built upon the ro- 
tational band in 163 Lu, which is identified with a rotation of a triaxial, strongly deformed shape. 
We find that the presence of pairing correlations substantially affects the energy of the wobbling 
excitations. Our calculations predict an onset of a tilted rotation at a critical rotational frequency 
where the energy of the wobbling excitations approaches zero. 
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The study of rotational and vibrational excitations is a 
major source of our understanding of a nuclear dynamics 
at low energy. For instance, the concept of deformation 
of a nuclear shape, created by an effective nuclear poten- 
tial, is well established. The analysis of interplay of a 
nuclear shape and orientation of an angular momentum 
led to various important discoveries [ij, |2j . Recently, ro- 
tating non-axially deformed nuclei attract a considerable 
experimental attention. For a long time, experimental 
data did not provide an irrefutable proof for existence 
of the non-axiality. The modern generation of detectors 
opens new avenues to study nuclear spectra and, in par- 
ticular, low- lying excitations near yrast line with high 
precision. In fact, the non-axial deformation gives rise 
to a new type of dynamics involving the orientation de- 
gree of freedom. This includes the chiral rotation |2| and 
wobbling excitations 0, 0, 0, H . 

As a rule, low-lying rotational bands are well described 
within the cranking model with a principal axis rota- 
tion. If a nuclear shape becomes non-axial, one expects 
low-lying vibrational excitations built on top of such a 
rotational band. These excitations, called the wobbling 
excitations by analogy with a classical wobbling motion, 
are created by a fluctuation of the angular momentum 
direction around one of the principal axes of a deformed 
nucleus. Having in mind that nuclei can rotate also about 
a tilted axis [2j, one may suggest for non-axial shapes a 
possible transition from a wobbling mode with a rota- 
tional axis fluctuating around one principal axis towards 
a tilted rotation with a fixed rotational axis lying in a 
plane in between two principal axes [g- 

We recall that the wobbling excitations in rotating nu- 
clei were suggested in Ref. |l| and analysed first within the 
microscopic approach in Ref. 0. However, the first ex- 



perimental evidence of such excitations is reported only 
recently i I M M M El 

The interest in wobbling excitations was sparked by the 
discovery of a particular excited rotational band above 
the rotational band in 163 Lu which is associated with 
a triaxial, strongly deformed (TSD) nuclear shape. The 
lowest band is called TSD1 band, while the excited one 
is denoted as TSD2 band. In accordance with the rates 
of the observed inter-band E2 transitions the structure 
of the TSD2 band is identified with wobbling excitations. 
Later, in the same nucleus a second rotational band was 
found, which was interpreted as a two phonon wobbling 
band Q. The fact that the two phonon band has less 
then twice the excitation energy of the one phonon band 
has been seen as a sign of non-harmonic vibrations [8llll|. 
The theoretical description of those band structures was 
done in terms of the phcnomenological, particle plus rotor 
model At the same time, a non self-consistent, 

microscopic analysis (within a mean field plus random 
phase approximation approach) were performed in Refs. 
ha and hj In Ref. hj it was concluded that the pair- 
ing correlations do not affect the wobbling phonon, and 
this should be considered as a specific feature of such 
excitations. On the other hand, it was also found that 
the wobbling motion is very sensitive to a single-particle 
alignment. It is well known, however, that the alignment 
decreases the pairing correlations. Therefore, the ques- 
tion arises about the validity of this conclusion. One of 
our goals is to clarify this issue in a calculation based on 
the cranking+random phase approximation (called here- 
after CRPA), with and without the pairing interaction 
and where we take care of the self-consistency. Following 
the analysis of Ref. ||| we also aim to study a possi- 
ble transition to a tilted mean field solution at the point 
where the RPA solution goes to a zero energy. 

We start with a pairing+QQ Hamiltonian [l8|: 
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Here are the single-particle energies of the spherical 
modified oscillator Hamiltonian h sp h 
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The operators c\ (c&) are fermion creation (annihila- 
tion) operators with the suffix k(l) labelling a com- 
plete set of quantum numbers. The parameters k and 
jS are standard ones [l^. The quadrupole residual in- 
teraction Q ■ Q = J2m=~2 Qm ^ S a Sum 0ver tne nve 

components of the quadrupole operators built up from 
f 2 Y 2 o and the linear combinations f 2 (Y 2m ± Y 2 _ m ) for 
m = ±(1,2). The quadrupole operators are defined as 
Qm = J2kir1m,kic{ T cir, (m = 0, ±1,±2) where q mi ki = 
(k\Q m \l) and t = ±1 distinguishes neutrons and pro- 
tons, respectively. The pairing operator P has a usual 
form = J2k>o ^kr^k • ^ ne m d ex & refers to the time 
conjugated state. 

To study rotational properties of the system we per- 
form the Legendre transformation into the rotating frame 



H 1 = H -uj- J 



(3) 



where J = (^J x , J y , J Z J is the angular momentum oper- 
ator and uj = u> (sinf? cos tp, sin (9 sint/j, cos 9) is the rota- 
tional angular frequency vector. Here, 10 is the magnitude 
of the rotational frequency and 9 and tp are Euler (tilt) 
angles of the cranking direction. It should be pointed 
out that we are not restricted to the signature symmetry 
used in the approach of Refs. andll7l 

The mean field part of if 7 , Eq. (|3J), can be written as 
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where A r are the pair field strengths and A T are con- 
straints for the particle numbers. Further, e 2 and 
7 denote the deformation parameters in the intrinsic 
frame of reference defined by the self-consistency condi- 
tions: k/Qq) = ^huj a e 2 cosj, n(Q 2 \ — ~^Tiuj e 2 sin7, 



Q±i,-2y = (see details in Ref. 0). Here < ... > means 
the mean field value. Using the Bogoliubov transforma- 
tion such as a\ = ^2 k Uikc\, + VikCk, we obtain Hartree- 
Bogoliubov equations that are solved for the TSD1 band 
in 163 Lu, with and without the pairing field. The major 
shells N — 4 — 6 are considered for both protons and 
neutrons. We used AiV — quadrupole interaction only. 
This type of interaction provides quite reasonable results 
with regard to equilibrium deformations in the moan field 
calculations (see also discussion in Ref. l2p|) . Further, this 
configuration space is sufficient for studying the pairing 
effects on low-lying excitations. The pairing strength is 
adjusted to reproduce the global fit of the odd-even mass 



differences for the ground state of even-even nuclei [2l| . 
The proton pair field is reduced by 20% due to the odd 
proton. The quadrupole interaction strength is fitted to 
reproduce the ground state deformation obtained with 
Nilsson-Strutinsky calculations. This is done separately 
with and without pairing. The self-consistent solution 
<f>" = |0), shortly denoted as ), corresponds to the min- 
imum of the energy surface E'(e 2l 7, A T , us, 9, <p) = (if 7 ). 
This implies that the equilibrium values for all these pa- 
rameters change as a function of the rotational frequency 
lu. Note, that our vacuum states ) are rotating odd- A 
particle configurations, similar to the ones used in Refs. 

Gain 

For huj < 0.5 MeV we find a principal axis rotation 
as the lowest solution. At Tilu w 0.25 (0.4) MeV the 
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FIG. 1: The neutron and proton pair field in the TSD1 band 
in 163 Lu as a function of the rotational frequency lu. 

proton (neutron) pair field disappears due to the grad- 
ual breaking of quasi-particle pairs (see FigHJ. We also 
find a increasing 7-deformation and a slowly changing 
62-deformation along the band, see Fig(2J At small rota- 
tional frequencies strong proton and neutron pair fields 
affect the deformation. In contrast to the unpaired case, 
the calculations with the pairing forces predict a small 
7-deformation at low rotational frequencies. Since the 
triaxial minimum is quite shallow in the unpaired calcu- 
lations, the presence of the pair field is enough to almost 
restore the axial symmetry. With the increase of the ro- 
tational frequency the equilibrium deformations manifest 
a triaxially, strongly deformed shape of 163 Lu. 

Once the self-consistent mean cranking solutions are 
found, we app ly the quasi-boson approximation in stan- 
dard way |18| in order to construct the vibrational wob- 
bling excitations by the RPA approach. The particulari- 
ties of this method for the rotational case can be found, 
e.g., in Ref. \^ The CRPA Hamiltonian is diagonalised 
by solving the equations of motion. As a result, we ob- 
tain the determinant of the secular equations (see details 
in Ref. |23|) which is solved numerically. 

The RPA equations have several spurious solutions re- 
lated to the symmetries, i.e., the rotational invariance 
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FIG. 2: The deformation parameters £2 and 7, associated with 
the TSD1 band in 163 Lu, as a function of rotational frequency 
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and the particle number conservation, broken in the mean 
field calculations. If the mean field problem is solved 
with a high accuracy, the spurious solutions connected 
with operators J 2 and N T will appear at zero energy and 
the solution connected with the operator T + ~ J x — iJ y 
will appear at the rotational frequency. Thus, the spuri- 
ous solutions are completely decoupled from the physical 
solutions. In contrast to the approach of Ref. H3, we 
obtain the RPA solutions related to the wobbling exci- 
tations from the full RPA determinant. In Ref. [17] this 
determinant is reduced to a simple dispersion equation 
for the wobbling excitations, which is valid , if and only 
if all spurious solutions are separated from the physical 
solutions. Accordingly, the numerical analysis of this dis- 
persion equation alone is not a warrant for a decoupling 
of the physical wobbling excitations from the spurious 
solutions. Due to the admixture of the spurious modes 
such an analysis would not be fully reliable. 

In Fig.[5]we compare our results of the wobbling excita- 
tions with experimental data Our results without 
the pai ring correlations are similar to the ones obtained 
in Ref. 16. The calculations predict an almost constant 
wobbling excitation energy up to Tiw « 0.4 MeV. Above 
this value we obtain a rapid decrease of the wobbling ex- 
citations, which leads to transition into a stable tilted 
solution. In contrast to the conclusion of Ref. 0, the 
pairing interaction dramatically changes the results at 
small rotational frequencies. The wobbling excitation 
energy is substantially larger at small rotational frequen- 
cies, while it is decreasing with the increase of the rota- 
tional frequency. The rate of the reduction as a function 
of the rotational frequency is slightly faster than it is 
seen in the experiment. After the collapse of both pair 
fields at hco « 0.4 MeV, the paired and unpaired calcu- 
lations predict similar results. Note, that for zero pair 
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FIG. 3: The wobbling excitation energy and tilt angle tp(8 = 
0) in the TSD1 band in 163 Lu as a function of rotational 
frequency uo. Results with and without pairing correlations 
are compared to the experimental data taken from j^Q. 



gap even though pairing vibrations are still present in 
RPA they do no longer mix with the quadrupole vibra- 
tions. Large values of the wobbling excitation energy at 
small rotational frequencies are brought about by strong 
pairing fields that reduce 7-deformation. These results 
confirm the prediction with regard to the 7 vibra- 
tional mode of the negative signature. According to Ref. 

the increase (decrease) of the rotational frequency (the 
pair field) transforms the negative signature 7 vibrational 
states with odd spins to the low-lying wobbling excita- 
tions with the increase of the triaxiality. 

It is a well known problem of cranked mean field cal- 
culations that the collapse of the static pair field at high 
rotational frequency happens too abrupt. This fact may 
explain why the reduction of the wobbling excitation en- 
ergy comes too fast in our calculations compared to the 
experimental data. As a consequence, we also obtain a 
transition to a tilted solution at a relatively low rotational 
frequency. This discrepancy could be resolved probably 
by introducing a self-consistent treatment with particle 
number projection. 

The ratio of the inter-band to intra-band E2 transi- 
tions for the lower part of the wobbling band is shown 
in Fig. 21 We use an effective charge of 0.5 (1.5) for 
neutrons (protons) to compensate the restricted size of 
our configuration space. One could observe a reasonable 
agreement with the experimental data, when the pair- 
ing correlations are included into the calculations, even 
though the calculated ratio is too large. Without the 
pairing correlations, we obtain a ratio which is smaller 
than is seen in the experiment. For comparison, we in- 
cluded the results from Ref. Il6| which are also smaller 
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FIG. 4: The ratio of the transition probability from the wob- 
bling band to the lower band and the intra-band transitions. 
The results with and without pairing are compared with the 
experimental data taken from 0, and with the results of 
Ref. Hrl We connect experimental points by a solid line for 
the sake of a comparison. 



than the experimental data. 

Summarising, the self-consistent treatment of the wob- 
bling excitations in the CRPA provides a good descrip- 
tion of the experimental data for 163 Lu. We conclude 
that the pairing interactions change the energy of the 
wobbling phonon, at least, indirectly by changing the 
7-deformation. In addition, our calculations indicate 
the onset of a tilted rotation in 163 Lu above a critical 
rotational frequency. The disappearance of the energy 
splitting between the two signature partner bands is one 
possible indication for a transition to a tilted rotational 
regime. Another indication would be that the two bands 
are connected via mixed Ml and E2 Al = 1 transitions. 
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